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Introduction: The dynamics of lava flow move- 
ment are controlled by the fluid interior. Crust, solids, 
and nondeformable material can only retard the ad- 
vance or spreading of a lava flow. Figure 1 shows a 
typical large, channelized lava flow found on the Mars 
plains. It has been suggested in [1] that such large lev- 
eed flows on the Mars plains were emplaced by a bal- 
ance between the formation and shedding of crust as 
the flow advances. For the prototypical flow north of 
Pavonis Mons (Fig. 1), such a balance leads to a flow 
morphology that approximately self-replicates at all 
locations along the flow path [2,3], Moreover, most 
quantitative characteristics of emplacement (e.g., vis- 
cosity, volumetric flow rate) of the prototype flow at 
Pavonis Mons resembled those of large channelized 
lava flows on Earth. The exception was the relatively 
long, sustained supply of lava, on the order of a year as 
opposed to hours or days for terrestrial analogs. 



Fig 1. A LEVEED LAVA FLOW NORTH OF 
PAVONIS MONS. 

Here we analyze and compare the quantitative em- 
placement characteristics of 5 additional large leveed 
lava flows north of Elysium Planitia, the plains north 
of Pavonis Mons, and near Arsia Mons. The issues 
investigated here are: 1) Are there significant regional 
differences in the viscosities, eruption durations, and 
volumetric flow rates of the large leveed lava flows? 
2) Are the emplacement mechanisms reasonable ex- 


trapolations of terrestrial experience? 3) Can averaging 
dimensional data provide an efficient method for re- 
gional comparisons of emplacement conditions? 

Approach: Figure 2 is a cartoon of the flow 

model in [3]. This model assumes steady-state up- 
stream conditions in the proximal zone and time- 
dependent levee-building only in the distal zone [2,3]. 
The essential inputs to the model are the channel width 
(w c ), levee width (w/), flow thickness (hi), the length of 
the flow (L), and the underlying slope (6). The outputs 
of the model are the thickness of the fluid interior (h L ), 
the thickness ( h cn , s , ) of the “crust” (interpreted as all 
nondeformable lava overriding the core), the time ( T ) 
required to form the crust (interpreted as the emplace- 
ment time), the average flow velocity, the viscosity (p ) 
of the molten core, and the volumetric flow rate (Q). 
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Fig 2. LEVEED FLOW CARTOON (FROM [3]) 


To continuously construct triangular shaped levees 
in the distal zone, the upstream core thickness must be 
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This criterion provides an immediate test of the valid- 
ity of steady-state upstream flow conditions for actual 
flows. Core thickness cannot be negative, nor can it 
exceed the overall thickness of the flow. These re- 
quirements give the validity condition 


0.1 9w <w t < 0.5w (2) 

When either of these limits is grossly violated for a 
particular flow, nonsteady upstream conditions must 
have prevailed to produce the observed deposit dimen- 
sions. Too much volume in the levee (according to eq 
2) indicates a time-dependent transient upstream con- 
dition such as a major overspill, a breakout, or a chan- 
nel surge due to either a significant upstream collapse 
or a spike in the effusion rate. 

Once the core thickness is known, the time required 
to construct the crust can be determined by the empiri- 



cal data analysis of [4] or other theoretical methods 
that account for exposures of the hot core due to crack- 
ing of the flow surface, stretching, shearing at the 
channel margins, etc. Using results from [4] for 
crustal thickening by conduction alone, the flow vis- 
cosity is estimated by 
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where k and X are empirical constants given in [4], 

Applications: Five additional long leveed flows 
with continuous channels over the length of the flow 
were identified in three, low slope plains settings. All 
flows have been studied in previous investigations 
[2, 3, 5, 6]. However, viscosity and flow rate estimates 
of past studies have been constrained only by 3 - 4 
orders of magnitude or more. 

The dimensions of the five flows all satisfy the va- 
lidity requirement in (2). It is noteworthy that several 
dozen candidate channelized flows [7] were investi- 
gated on the steeper slopes of Ascraeus Mons. How- 
ever, none satisfied the validity requirement shown in 
(2) suggesting that time-dependent effects induced by 
variable eruptions rates or topographic slope changes 
were significant factors. 

Baloga and Glaze [2,3] applied the levee model to 
multiple topographic profiles from MOLA gridded 
data for the prototype PavE flow (Fig. 1), then per- 
formed a somewhat complex statistical analysis to ex- 
tract the systematic behavior of the levee construction 
and flow dynamics. Here, the procedure is reversed, 
i.e., average dimensions were obtained first and the 
levee model was applied subsequently. The averages 
used in the computations for PavE and the five addi- 
tional flows are shown in Table 1 . 

Table 1. 


Flow 

h m 

w r km 

w/ km 

L km 

6 deg 

PavE 

47 

14 

12 

173 

0.05 

PavW 

32 

23 

9.0 

307 

0.09 

Ely 13 

66 

2.7 

1.0 

43 

1.08 

Ely 12 

54 

3.7 

1.2 

107 

1.08 

Ely 9 

58 

5.3 

1.4 

107 

1.08 

Arsia 

22 

1.9 

0.8 

75 

0.7 


Quantitative inferences for the emplacement condi- 
tions are shown in Table 2. Processes such as over- 
spills, surficial cracking, internal flow circulations and 
shearing at the margins, and responses to small-scale 
topography typically change results by factors of 2 to 
4. Therefore, the results for Table 2 have been rounded 
to reflect these uncertainties in the actual emplacement 
conditions. Results for PavE were taken from the pro- 
file-by-profile analyses in [2,3] for comparison. 


Table 2. 


Flow 

T d 

O m3/s 

p Pa-s 

PavE* 

600 

10 3 

10 6 

PavW 

700 

10 3 

10 4 

Ely 13 

2700 

10' 

10 7 

Ely 12 

1000 

10 2 

10 6 

Ely9 

300 

10 3 

10 6 

Arsia 

400 

10 2 

10 s 


Although flow lengths and thicknesses vary by fac- 
tors of 3-4, slopes a factor of 20, and widths a factor of 
1 0, there are a number of off-setting factors in the dy- 
namics of this model. This was noted in [3] as a poten- 
tial explanation why channelized lava flows over many 
scales of dimensions appear similar on different plan- 
ets and in dramatically different volcanologic settings. 

The long channelized plains flows have somewhat 
remarkable similarities in the emplacement dynamics. 
The volumetric flow rates are at the high end range of 
terrestrial experience for large basaltic eruptions. The 
viscosities of the fluid cores similarly are comparable 
to terrestrial estimates of basalt flows in distal seg- 
ments. The notable difference, however, is the duration 
of the supply, persisting for approximately 1-10 yrs as 
opposed to hours or days on terrestrial shields. 

Conclusions: This levee flow model provides a 
means for isolating the viscosity and volumetric flow 
rate instead of the product of the two. Emplacement 
times can be significantly longer than results by other 
models because the dynamics are controlled only by 
the inner fluid core, which may be significantly less 
than the total flow thickness. It appears that cautiously 
applying the model to averaged dimensional data is an 
efficient method for establishing regional differences 
or similarities of emplacement conditions. 

Large channelized flows on the Mars plains have 
relatively high volume flow rates and core viscosities 
by terrestrial standard, but the values are not at all un- 
reasonable. The key difference in all regions investi- 
gate to date is the duration of supply. This has signifi- 
cant implications for subsurface conditions and proc- 
esses that can sustain such conditions. 
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